Testing General Relativity with Current Cosmological Data 
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Deviations from general relativity, such as could be responsible for the cosmic acceleration, would 
influence the growth of large scale structure and the deflection of light by that structure. We clarify 
the relations between several different model independent approaches to deviations from general 
relativity appearing in the literature, devising a translation table. We examine current constraints 
on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, 
cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. 
A Markov Chain Monte Carlo likelihood analysis of the parameters over various redshift ranges 
yields consistency with general relativity at the 95% confidence level. 
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I. INTRODUCTION 

The nature of gravitation across cosmological ages and 
distances remains a frontier of current knowledge as we 
try to understand the origin of the cosmic acceleration 
1 , |2[ . Newly refined observations of cosmic structure 
3, |j] make it possible to test the predictions of general 
relativity (GR) for its influence on the growth of cos- 
mic structure through gravitational instability and the 
gravitational lensing deflection of light by that structure. 
Indications of a deviation from GR would have profound 
consequences for cosmology, as well as for fundamental 
physics. 

To explore for new gravitational phenomena, it is use- 
ful to parameterize the deviations from GR in the grav- 
itational field equations. A common approach is to in- 
troduce two new parameters. The first parameter im- 
poses a relation between the two gravitational potentials 
entering Newton's gravitational law of acceleration and 
the Poisson equation. These are equal in GR in the ab- 
sence of anisotropic stress but different in many theories 
of modified gravity. The second parameter establishes a 
new relation between the metric and matter through a 
modified Poisson-Newton equation, which can be viewed 
as turning Newton's gravitational constant into an effec- 
tive function of time and space. Numerous realizations 
of these relations have been put forward in the literature 

One motivation for our study is to attempt to relate 
these disparate, but closely related, approaches. Further- 
more, many studies have focused on the ability of future 
measurements to discriminate among various models and 
to carry out parameter estimation [19l - l3Cf , however there 
is sufficient data at present to evaluate preliminary tests 
of GR [3ll - l34l . [3fj| - l38j . We concentrate here on current 
constraints, which also allows us to examine a recent 
claim of a possible departure from GR \^ . 

The main points of this article are thus to 1) clar- 



ify the relation between different parametcrizations and 
what the degrees of freedom are in a consistent system of 
equations of motion, 2) confront the parameters encoding 
deviations from GR with current data to test the theory 
of gravity, and 3) discuss which features of the data have 
the most sensitivity to such a test and what astrophysical 
systematics may most easily mimic a deviation. 

In Sec. we lay out the gravitational field equations 
in terms of the metric potentials and matter perturba- 
tions and compare several forms of parameterizations, 
giving a "translation table" between them. We illustrate 
in Sec. IHII the influence of the parameters on the cosmic 
microwave background (CMB) temperature power spec- 
trum, the matter growth and power spectrum, and the 
weak lensing shear statistics. Using Markov Chain Monte 
Carlo (MCMC) techniques, we then constrain the devia- 
tion parameters with current data in Sec. IIVI We briefly 
discuss astrophysical systematics and future prospects in 
SecEl 



II. SYSTEMS OF PARAMETERIZING 
GRAVITY 

The most accurate observations of the effects of grav- 
ity have been made in the local universe, e.g. within the 
solar system and in binary neutron star systems fidl - |43T ] . 
These observations can be used to distinguish between 
various theories of gravity throu gh the parameterized 
post-Newtonian (PPN) formalism pdl |45| . The standard 
PPN formalism introduces a set of constant parameters 
that take on various values in different gravity theories. 
This, however, does not give a full description of possi- 
ble deviations from General Relativity over cosmological 
scales. 

Recent interest in modified gravity has concentrated 
on those theories that can serve as an alternative expla- 
nation for the current period of accelerated cosmic ex- 
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pansion. In order for modifications producing late-time 
acceleration on cosmic scales to agree with local tests 
of gravity they must contain length and/or time depen- 
dent modifications, which do not occur in the standard 
PPN formalism. Moreover, for some theories the nat- 
ural arena for the PPN formalism - solar system and 
binary neutron star system observations - may be less 
discriminating than cosmological tests of gravity, given 
that the modifications are on large scales. This has led 
to efforts to establish a parameterized formalism that al- 
lows for meaningful comparison between modified gravity 
theories within a cosmological framework [5l-fl8j , without 
assuming a specific model. 

A. Degrees of Freedom 

Changes in the laws of gravitation affect the relation- 
ship between the metric and matter variables. Let us 
explore the degrees of freedom available to define this 
relation. Restricting our attention to scalar degrees of 
freedom of the gravitational field, the metric has only 
two physically relevant scalar functions, or potentials, 
given by the line element (in conformal-Newtonian gauge, 
adopting the notation of jir| ) 

ds 2 = a 2 [-(1 + 2V>) dr 2 + (1 - 20) dx 2 } , (1) 

where a is the scale factor, r the conformal time, and x 
the spatial coordinate. In addition to the metric poten- 
tials 4> an d ijj, perturbations to a perfect fluid introduce 
four additional scalar functions: density perturbations 
dp, pressure perturbations Sp, velocity (divergence) per- 
turbations 0, and a possible nonzero anisotropic stress 
a. 

The dynamics of any particular theory are then spec- 
ified when six independent relations between these six 
quantities are given. Further restricting attention to 
those gravity theories that maintain the conservation of 
stress energy, V M T^„ = 0, the resulting generalized conti- 
nuity and Euler equations give two scalar equations and 
the grav itational field equations supply the remaining 
four gfij. 

Since the cosmic expansion shifted from deceleration 
to acceleration only recently, since z < 0.5 0, grav- 
ity theories that account for this transition without any 
physical dark energy require a significant departure from 
GR at late times. Consequently, nonrclativistic matter 
is the dominant component of the cosmological fluid and 
so 5p — Sp m = and a = a m = 0. Hence, in these 
theories the dynamically important equations consist of 
two, as yet unspecified gravitational field equations and 
the two equations of stress-energy conservation applied 
to matter, which in Fourier-space are given by 

L = -e m + 30, (2) 
e m = -w m + k 2 ip. (3) 

In the above equations, 5 m = Sp m /p m with p m the ho- 
mogeneous part of the matter density, H = a/a, the dot 



denotes a derivative with respect to conformal time, and 
k is the wavenumber. There still remains freedom in 
setting the two gravitational field equations to close the 
system, subject to the requirement that the theory ap- 
proaches GR within the solar system. 

The two field equations that can close the system in 
the case of GR are 

V 2 = 47rGa 2 p m A m , (4) 

ijj = <\>, (5) 

where 

A,„ = S m + -p-flm ■ (6) 

In a wide variety of alternative theories of gravitation, 
additional scalar degrees of freedom modify the strength 
of Newton's constant, and enforce a new relationship be- 
tween the potentials 4> and ip. Therefore, one choice for 
the modified field equations in Fourier-space is 

- k2 A A + + B B = ^Gp(r,k)p m A m , (7) 
<j> = rj(T,k)iP, (8) 

where A and B are constants, and p and rj are func- 
tions of time and scale, which are still to be determined. 
As we will see, there are many other choices that can 
be made for the exact form of parameterization. These 
choices influence the constraints and the correlations be- 
tween those constraints that particular observations give 
for a particular set of post-GR parameters. We discuss 
some of the frameworks in the next subsections. 



B. tc/iCDM 

We refer to the equations of motion used in 0, HH, 
H3, [37| as zuCDM. In tuCDM, the equations of motion 
for cosmic perturbations are determined by enforcing the 
relation 

i> = [1 + w(t, k)] 4> (9) 

for the potentials arising from non-rclativistic matter, 
where the departure from GR is controlled by the pa- 
rameter w. In practice, this is carried out by adding 
a source to the off-diagonal space-space Einstein equa- 
tion in order to simulate a smooth transition from GR to 
modified gravity. 

Next, requiring that the new gravitational phenom- 
ena do not introduce a preferred reference frame distin- 
guished by a momentum flow, e.g. a 9 that would be at- 
tributed to a dark fluid, the time-space Einstein equation 
is preserved, whereby 

- fc 2 (J> + Hip) = -4irGa 2 {p + p)6 . (10) 
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As discussed in [33], preserving the time-space Einstein 
equation along with the modification in Eq. ([3]) still re- 
sults in a correction to the GR Poisson equation. This 
can be thought of as a consequence of the conservation of 
stress-energy and the related Bianchi identity as applied 
to the modified gravitational field equations. 

We now propose to extend the roCDM model, to incor- 
porate a new parameter fi that controls the modification 
to the Poisson equation, i.e. 



tu/xCDM parameters as follows: 
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D. Gravitational Growth Index -yc 
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— k 2 <p = /i(r, k) AnGa 2 p m A T 



(11) 



Procedurally, this equation replaces Eq. (|10p for obtain- 
ing the evolution of the gravitational fields. We call this 
new parameterization zupCDM and note that in this pa- 
rameterization the time-space Einstein equation is gener- 
ally modified, as opposed to in the wCDM parameteriza- 
tion. Note that setting p = 1 does not reproduce the orig- 
inal znCDM model since there vj itself modifies the Pois- 
son equation as discussed above. This parametrization is 
consistent with the conservation of large-scale curvature 
perturbations following an argument made in Ref. [2l| . 
They argue that super-horizon curvature perturbations 
are conserved, so long as the velocity perturbation 9 is of 
order (k/H) 2 . This can be seen to be true from Eq. (|A3[) 
presented in the appendix of this work. 



C. PPF Linear Theory 

A parameterized post-Friedmann (PPF) framework of 
linear fluctuations was introduced by 0, [lj| to describe 
modified gravity models that yield cosmic acceleration 
without dark energy. It captures modifications of gravity 
on horizon, sub-horizon, and non-linear scales. Once the 
expansion history is fixed, the model is defined by three 
functions and one parameter, from which the dynamics 
are derived by conservation of energy and momentum and 
the Bianchi identities. Modifications to the relationship 
between the two metric perturbations are quantified by 
the metric ratio 
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In the linearized Newtonian regime, a second function 
/c(a) relates matter to metric perturbations via 



k 2 {cj) + 4>) 
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a 2 p r , 



(13) 



The corresponding quantity that defines this relationship 
on super horizon scales is f((a). The last quantity that 
needs to be defined is cr, which determines the transition 
scale from superhorizon to quasistatic behavior in the 
dynamical equations (see 0, [l3| for details) . 

The PPF parameters can be directly related to the 



Another way to close the system of equations is to spec- 
ify the evolution of one of the perturbed fluid or metric 
variables. A standard choice is to determine a specific 
evolution for A m through the gravitational growth index 
7G introduced to parameterize deviations from general 
relativity in growth by [Jj| • This was partly tied to the 
metric potentials in [49( but here we present a more com- 
plete relation. 

From Eq. (23) of [49| we see the key quantity is the 
modification of the source term in the Poisson equa- 
tion, there called Q. The second order equation for 
the evolution of the density perturbation arises from 
"V 2 ip, and there is also a modification p allowed in 
the gravitational coupling as in Eq. (jTTJ) . In essence, 



-fc 2 (l + w)(j) (1 + x 4nGa 2 p m A m . Thus 



Q = (1 + w)p. The relationship between w, p and the 
evolution of A m is presented rigorously here in Eq. (|A5|) 
(also see Sec. lIIIBI) . 

The gravitational growth index in Eq. (23) of [4£| thus 
relates to the vupCDM formalism through 
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Note u>oo is an effective high redshift equation of state 
defined in terms of how the matter density in units of the 
critical density, fi m (a), deviates from unity (specifically, 
Woo = [dlnJ7 m (a)/dlna]/[3(l — il m (a)]). In the last line 
of Eq. (|T7|) we specialize to a ACDM expansion history, 
as used throughout this article, so Woo = — 1, and to the 
ansatz for vo and p, used later in Eqs. (Il£ 



E. Relating Parameterizations 

The discussion above is by no means an exhaustive list 
of the parameterizations proposed in the literature to de- 
scribe departures from GR. Many more exist, and while 
all of them have in common a relatively simple param- 
eterization of the departure from cf> = ij), they all differ 
in how they close the system of equations. Some, like 
tu/iCDM, modify the Poisson equation directly. Others, 
like ijjCDM, retain one of the Einstein equations. 

Table U lists some of the most common parameteri- 
zations and presents a useful translation between their 
post-GR parameters and wpGDM. With the possible 
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exception of the parameterization from [39( (see next 
paragraph and its footnote), all of the parameterizations 
presented are presumed to leave the equations of stress- 
energy conservation unmodified. 

Since none of these model-independent approaches 
start from an action, one must be careful to trace the sys- 
tem of equations to make sure that the phenomenological 
modifications do not under- or over-specify the system 1 
and do satisfy stress-energy conservation. Another ap- 
proach involves testing consistency relations valid in GR 
between observables; see for example @, HH, H3|. 



III. INFLUENCE OF GRAVITY 
MODIFICATIONS ON OBSERVATIONS 

The behavior of the CMB, weak lensing, and matter 
power spectrum in the -zuCDM scenario have been dis- 
cussed in [H, H3] ■ The consequences are slightly different 
when we introduce fj. in the nj/iCDM parameterization. 
In the case that w < and fx < 1, both lead to an am- 
plification of low-^ CMB power; w > and /j, > 1 both 
suppress it. This allows us to play the two parameters 
against each other, combining positive (negative) values 
of vj with smaller (larger) values of /i to generate non- 
GR power spectra that appear to be in better agreement 
with the data than those obtained within the confines of 
the tuCDM model. That either parameter can enhance 
or suppress power results in a degeneracy between w and 
/i in any multi-parameter exploration of the data. Obser- 
vations that can break this degeneracy therefore become 
vital to diagnosing departures from GR. 

For the purposes of the discussion in this section, we 
will assume the redshift dependences 

3 

ZD = VDqCC 

/i = l + /i a 3 . (18) 

Note this form can be motivated by the scaling argu- 
ment in [49j , that the deviations in the expansion history 
should keep pace with the deviations in the growth his- 
tory. Otherwise one tends to either violate GR at early 
times (causing difficulties for primordial nucleosynthesis 
and the CMB) or does not achieve acceleration by the 
present. In addition to the CMB, we also discuss the 
effects of our post-GR parameters on the matter power 
spectrum and on weak lensing statistics. 



A careful reading of [39ll reveals that there four unknowns - <f>, 
ip, S, and 9 - arc evolved with five equations — the continuity 
equation Eulcr equation Poisson equation Q, and the 
post-GR parameter equations </> = r/ip and A m = MA ra fim . 
Thus, the system is over-specified. 



A. CMB Anisotropy Spectrum 

We modified versions of the public Boltzmann codes 
CMBfast [5l| and CAMB [HJ to evolve the cosmological 
perturbations according to parameterization (|18[) and the 
equations of motion presented in Sec. iBBl We used these 
codes to generate examples of CMB anisotropy and mat- 
ter power spectra for different values of zuq and /io; in 
order to focus on the non-GR effects, in this section all 
other cosmological parameters are set to their WMAP5 
maximum likelihood values [53j. Figure [J shows the re- 
sulting CMB anisotropy spectra. As in [34| for tuCDM, 
negative values and extreme positive values of the post- 
GR parameters amplify the power in the low-£ multi- 
poles. Moderate positive values suppress the low-£ power. 
This is a manifestation of the integrated Sachs- Wolfe 
(ISW) effect. The high-£ power is unaffected. 

The ISW effect arises when time evolving <p and ip po- 
tentials cause a net energy shift in CMB photons. The 
CMB ISW power is sourced as 

C,~(^ + ^) 2 . (19) 

As was discussed in [33], the evolution of <p an d ip po- 
tentials in the universe is a competition between gravita- 
tional collapse trying to deepen the potentials and cosmic 
expansion trying to dilute them. Under GR with a cos- 
mological constant, the expansion wins and the source 
term for the ISW cp + tp > (note cf>, ip < 0). By weaken- 
ing gravity, wq or /^o < tilts the competition even more 
towards cosmic expansion, hastening the dilution of <p 
and ip, causing <p + ip to be even larger, and amplifying 
the ISW effect. Positive wq or /i amplifies gravity - ei- 
ther by directly deepening the Newtonian potential ip so 
that mass is more attractive (wq > case) or by causing 
A m to source a deeper potential through the modified 
Poisson equation (/io > case) - so that the dilution due 
to cosmic expansion is slowed, leading to a weaker ISW 
effect. In the case of extremely positive zuq or the 
ISW deepening is so pronounced that the sign of (p + ip is 
reversed, but since the ISW effect in the power spectrum 
depends on the square, the ISW effect is again ampli- 
fied. High-^ power is unaffected because the ISW is a 
sub-dominant effect on those scales. 

Figure [2] more clearly illustrates this bimodal behavior 
by plotting the change in quadrupole power relative to 
GR as a function of the post-GR parameter, varying one 
at a time (compare Fig. 4 in [13] )■ The blue, dot-dashed 
curve is generated by varying zuq and holding fixed /io = 
0. The red, dashed curve is generated by varying ^iq and 
holding fixed vjq = 0. Note that the CMB appears to 
be more sensitive to differing values of /io than of wq. 
The black, solid curve is generated by varying wq and 
compensating for this by setting /io = 2/(2 + toq) — 1. 
This choice is motivated by the alternative definition of 
the unmodified Poisson equation 

- fc 2 (0 + = 47rGa 2 p m A m (20) 
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parameterization 



parameter 
relating cj> and ip 



closing 
parameter 



comments 



zuCDM [6, 34, 37] w 
Curvature [12] "/bz ~ 

PPF[9,PJ g 



2 + 137 



MGCAMB [20, 21] 7 mgc 

Cf. [7,8,i0] 7): 



i 

l+ro 



g 



Sub-horizon [101 
Growth index [48] 
Decoupled [39J 7) 



1 + VJ 
additional 



Retains equation (j 1 Of) 

Conserves curvature perturbations ( Effectively retains Eq. (flTJ 

See appendix in Ref. [34] 



n(2+t7j) 



1 



f^MGC = m(1 + ra?) 

fj _ _ f (2 + ro) 

r</> = m 
7g = it ( : 



1 + TJ7 



7G = 



ln(A m /«A m ) 
lnf2 m (a) 



Includes scale-dependent transition 
between super- and sub-horizon regimes 

Modifies Poisson equation with ip 
instead of <f> in Eq. (jl l|) 



Only defines (zu, /i) — > 7g not inverse 



Over-specified (also enforces Poisson eqn). 



TABLE I: Translation between several different parameterizations of modified gravity and the kj^iCDM framework. 



(see further discussion in the next section). We see that, 
for a wide range of values of zjq, complementary (vjq > 
and no < or vice-versa) values of /j,q cancel out much 
of the late-time ISW effect found in Fig. [TJ as alluded to 
in the introduction to this section. 



B. Matter Power Spectrum and Weak Lensing 
Statistics 

We investigate the power spectrum of the matter per- 
turbations 5 m as a function of wavenumber k in Fig.[3]for 
the same set of models. Again the most dramatic post- 
GR effects occur at large scales. This is not due to any 
scale dependence in the modifications (we took w and /i 
to be independent of k), but simply from the k 2 factor 
in the modified Poisson equation (jlll) . 

For the weak lensing shear correlation function, as for 
many other observables, we need to know how overdensi- 
ties grow with scale factor. In the case of GR and roCDM, 
this is a relatively simple proposition since the growth of 
overdensitics A m is scale-independent after decoupling. 
As just discussed, this no longer holds for ro/xCDM. It is 
possible, using energy conservation and Eqs. (|9|) and (|TTj) . 
to derive a second-order differential equation for the evo- 
lution of A TO . We show the derivation and result in the 
Appendix, and focus here on the parameter dependence. 

With the exception of one term on the middle line of 
Eq. (|A5j) . all of the terms containing metric potential 
modifications to general relativity (the fi and w terms) 
are multiplied by a factor of TL 2 /k 2 . Hence, we expect 
that the strongest departures from GR predictions oc- 
cur for small values of k. Since the most important as- 
pect for comparing modifications against observations is 



the change in shape of the power spectrum, rather than 
its normalization, in Figure |3] we normalize the power 
spectrum to agree with Pgr at large k. The strongest 
deviation in shape indeed occurs for k < 0.002 Mpc -1 . 
The one exceptional term in Eq. (|A5[) is precisely the 
(1 +w)fx term discussed in Sec. Ill Dl entering the gravita- 
tional growth index jg formalism, and this will dominate 
for large values of k, giving a scale- independent enhance- 
ment (suppression) for positive (negative) -cuq or /in- 

Figure |4] plots £_e, the E mode of the weak lensing 
shear two-point correlation function (Eq. (8) of Ref. Q), 
normalized to the value under GR as a function of an- 
gular separation on the sky. For the angular scales of 
interest, the effects of changing /Jo and zuq principally 
manifest themselves as a renormalization of £e. This 
is because the scales plotted are much smaller than the 
scales (k <~ %) at which shape-changing effects mani- 
fested themselves in Figure[3] Non-linear power is treated 
using the usual subroutine halofit based on the semi- 
analytic fitting scheme presented in Ref. (HJ. While we 
acknowledge that this is not strictly appropriate for mod- 
ified gravity, we have no reason to think that the effect 
will be substantial for reasonable values of /_io and wq. 
Furthermore, the constraints presented below in Section 
IIVI appear to principally derive from effects at the low-fc, 
rather than the high-fc, limit. 



IV. CONSTRAINTS ON DEVIATIONS FROM 
GR 

We now examine constraints imposed by current data 
on deviations from GR, allowing a large set of cosmo- 
logical parameters to vary simultaneously. The investi- 
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FIG. 1: CMB anisotropy spectra are plotted as a function 
of the parameters wo and /io in Eqs. (|18[) . As in [33], the 
post-GR effects all occur in the low-i? multipoles. The CMB 
anisotropy is more sensitive to variations in /to than to varia- 
tions in zuq. See Fig. [2] for more on this point and on varying 
wo and /io simultaneously. 



gation includes two different functional dependences for 
the gravitational modification parameters w{a) and /u(o). 
The first model for the post-GR parameter form does not 
assume a particular rcdshift dependence but allows w 
and p, to take independent values in each of three rcdshift 
bins. (In fact, we slightly smooth the transitions so as to 
avoid infinities in the derivatives entering the ISW effect, 
with a transition modeled by an arctan form of width 
Aa = 0.01.) That is, fi = {1 + fJ>oai 1 + £*06j 1 + Moc} an d 
zj = {■cuo a ,zuob,'cuoc} for {2 < z < 9, 1 < z < 2, z < 1}. 
We assume that w and [i are scale-independent. For 
z > 9 we assume that differences from GR are negligible 
so /i = 1 and zu = 0. 

We test this theory against the data using a mod- 
ified version of the public MCMC code GOSMOMC 
with a module (first presented in [5fj| ) to in- 
corporate the COSMOS weak lensing tomography data 
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® = 0, n = x 
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x (post GR parameter) 

FIG. 2: The change in quadrupole power relative to the 
value in GR is plotted as a function of roo and /io. The blue, 
dot-dashed curve shows the effects of varying wo with fixed 
/io = 0. The red, dashed curve shows the effects of vary- 
ing /io with fixed wo = 0. One can mimic the unmodified 
GR CMB spectrum over a much wider range of post-GR pa- 
rameter values by simultaneously varying wo and /to in op- 
posite directions, as shown in the black, solid curve using 
/to = 2/(2 + wo) — 1. The horizontal dotted line denotes 
perfect agreement with GR. 



and data from the CFHTLS survey 0. We also in- 
clude WMAP5 CMB data [57l459j and Union2 supernova 
distance data [6(J. In all cases, we use the full covariance 
matrix (including systematics in the Union2 case) pro- 
vided by the group who collected and initially analyzed 
the data. In addition to the post-GR parameters, the 
parameter set includes r2&/i 2 , Qch 2 , 9 (the ratio of the 
sound horizon to the angular diameter distance to last 
scattering), r (the optical depth to rcionization), n s , the 
amplitude of the SZ effect, and the amplitude of primor- 
dial scalar perturbations. We assume that w = — 1 for 
our effective dark energy, that = 0, and that there are 
no massive neutrinos contributing to dark matter. Each 
weak lensing data set requires 3 nuisance parameters. 
Thus, we integrate over a total of up to 16 parameters, de- 
pending on the data sets used and the parametrization of 
/i and w chosen. Under the binned parametrization, we 
vary /x or w but not both simultaneously, which would re- 
quire 19 independent parameters. This choice was made 
both for convenience and to reproduce the analysis of 
Ref. [HI. Additional MCMC calculations done in which 
both /x and w were allowed to vary in all three bins also 
returned results consistent with GR in the presence of a 
cosmological constant. 

Figure [S] shows the marginalized probabilities on the 
tt0a.b,c parameters for runs in which /ioo,&.c = 0, so that 
/i = 1 and the Poisson equation defined as in Eq. ([J| 
remains valid at all redshifts. Figure [6] shows similar 
constraints on fio a ,b.c in the case where vao a ,b.c = . Our 
results in all cases are consistent with GR within the 95% 
confidence limit, although they do allow the possibility 
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FIG. 3: We plot the matter power spectrum (normalized to 
k — IMpc" 1 ) generated by varying the parameters zuo and 
fio- Unlike under zuCDM [34|, even our scale-independent 
parameterization has scale-dependent effects due to the k 2 
factor in the Poisson equation. The bottom panel shows the 
residuals of the top panel, i.e. the deviation relative to GR 
when varying zuo (the fio case looks similar), to highlight the 
scale-dependent regime at low k and scale- independent regime 
at high k. 



of departures from GR with wq or /iq ~ 0.1. 
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FIG. 4: We plot the ratio of the E mode of the weak lensing 
shear two-point correlation function (Eq. 8 of H|) to the same 
statistic calculated in GR, with all parameters but either zuo 
or fio fixed, to see the influence of the non-GR parameters. 
For the most part, post-GR parameters serve to renormalize 
the correlation function. As with the CMB anisotropy and 
matter power spectra, the effect is more sensitive to changes 
in no than to changes in zuo- 



Constraints on the usual cosmological parameters are 
largely unaffected by the introduction of fi and zu. Mean 
values shift by less than ~ la and marginalized uncer- 
tainties are comparable between GR and non-GR MCMC 
runs. The only notable exceptions are as and fl c h 2 (the 
physical density of cold dark matter in the universe), 
whose marginalized uncertainties increase by up to a fac- 
tor 2.3 upon the introduction of post-GR parameters. 
This is consistent with the observation that fi and vu prin- 
cipally modify the growth history of cosmological pertur- 
bations. 

Figure [7J plots the 2-dimensional confidence contours 
for the post-GR parameters vuq, fio m the case of rcdshift 
dependence as in Eqs. (fT5|) . Note that since this param- 
eterization has the strongest effect at low redshift, the 
greater sky area of CFHTLS has more leverage in con- 
straining the parameters than the greater depth of COS- 
MOS. For the binned parametrization, the constraints 
from MCMC runs with WMAP5+Union2+CFHTLS (no 
COSMOS) were indistinguishable from those including 
COSMOS as well, supporting the supposition that the 
sky coverage of CFHTLS is, for current data, more im- 
portant than the redshift depth of COSMOS in constrain- 
ing the post-GR parameters. 

Table [IT] presents the 95% constraints on our post-GR 
parameters for all of the MCMC calculations considered 
in Figures [5][7l All of the results are consistent with GR. 

We also note that in Figure [7] the contours exhibit the 
same degeneracy implied by Figure [2] Apparently, the 
probe of growth provided by current weak lensing data 
is not able to add much more leverage to the CMB data. 
This can also be seen in the lack of significant change 
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Binned -ca, fj, = 1: Binned /i, zu — 0: Parameterization {TSjl: 

COSMOS +CFHTLS +CFHTLS COSMOS +CFHTLS 

-0.11 < zuoa < 0.12 -0.15 < vj 0a < 0.060 -0.074 < fi 0a < 0.080 -1.4 < ro < 2.8 -1.6 < ro < 2.7 

-0.098 < ro ( ,6 < 0.23 -0.13 < zu ob < 0.18 -0.058 < nob < 0.14 -0.67 < /i < 2.0 -0.83 < fia < 2.1 

-0.054 < vjqc < 0.39 -0.074 < txj 0c < 0.33 -0.023 < jioc < 0.22 

TABLE II: 95% confidence limits on post-GR parameters in the MCMC calculations considered in Figures \E\ (left two columns), 
H (middle column), and (right two columns). Columns labeled "COSMOS" use WMAP5+Union2+COSMOS data. Columns 
labeled "+CFHTLS" use WMAP5+Union2+COSMOS+CFHTLS data. Recall that in the binned parameterization, redshift 
bin a is 2 < z < 9, redshift bin b is 1 < z < 2, and redshift bin c is z < 1. 



in the width of the probability distributions in Figure [5] 
when adding weak lensing. 

The degeneracy illustrated in Figure [5] is plotted as the 
black, solid curve in Figure [7J The agreement with the 
likelihood contours is quite interesting, calling to mind 
the discussion in Sec. IIII A"l about parameter covariances. 
This arose from the observation that an unmodified Pois- 
son equation ([201) that relates the sum of the two metric 
potentials to the underlying density fluctuations leaves 
the large-scale CMB predictions nearly unchanged when 
varying the ratio of the metric potentials, i.e. w. That 
degeneracy is due to the fact that the large-scale CMB 
predictions depend on the sum of the two metric poten- 
tials (cf. Eq. [T9|) . If this sum is directly related to the 
underlying density perturbation then the only effect w 
can have on the large-scale CMB is through its effect on 
the evolution of A m ; by contrast, if the Poisson equation 
is of the form of Eq. ([TTj) . where only one potential is 
related to A m , then zu also appears in a multiplicative 
factor. Thus the specific approaches to modifying grav- 
ity give distinct relations between the parameters and the 
observables. 

For observations that depend on the combination <p+ij) 
there will be a degeneracy along the curve (see Eqs. fLU 
and [l"5j) 



We find numerically that this degeneracy applies approx- 
imately to both large-scale CMB as well as weak lensing 
observations, even though both measurements have a fur- 
ther dependence on zu and /i through the growth factor 
(cf. Eq. [17]). The relation in Eq. ([2TJ) gives the black, 
solid curve in Figure [7] and indeed is quite close to a 
degeneracy in the constraints. 

V. DISCUSSION 

Testing general relativity on cosmological length scales 
is an exciting prospect enabled by improvements in data. 
To interpret such a test requires an approach to pa- 
rameterizing modifications from GR, similar to the PPN 
method for tests within the solar system and using binary 
pulsars, but appropriate for cosmic scales. Numerous 
paramctcrizations have been suggested and we compare, 



and in some cases, unify them through a "translation" 
table. These approaches can effectively be interpreted 
within one formalism with two parameters zu and /i (an 
extension of the previous roCDM scenario). 

In this generalized tn^CDM model, even if the two pa- 
rameters characterizing modifications to gravity are scale 
independent we find effects that are visible in the large- 
scale structure matter power spectrum, and thus in weak 
lensing shear correlations, that depend on scale. We give 
quantitative results for the effects of the modifications 
on the cosmic microwave background temperature power 
spectrum, the growth of matter density perturbations 
and the density power spectrum, and the weak lensing 
statistics, along with analysis of the physical basis of the 
effects. On large scales in the density power spectrum, 
values of zu or /x above their GR values cause suppression 
of power while leading to enhancement on smaller scales. 

We confront the modifications to GR with current cos- 
mological observations, analyzing CMB (WMAP 5-year), 
supcrnovac (Union2), and weak lensing (CFHTLS and 
COSMOS) data. Employing two different forms of de- 
pendence of the modifications on redshift, we find no ev- 
idence at 95% confidence level for such extensions to GR, 
regardless of the combinations of data used. Note that 
this holds for both the data employed by [39| (which used 
an overspecified system of equations in that analysis), 
and a more comprehensive set of observations. 

We also verify the trade-off between zu and fj, pre- 
dicted analytically. Such covariance leads to an inter- 
esting degeneracy for measurements depending on the 
sum of the metric potentials, although growth measure- 
ments depend on a different combination. Since large 
scale CMB and weak lensing depend on the sum of the 
metric potentials, one could consider the Poisson equa- 
tion for the sum, and here the key parameter is the ef- 
fective Newtons constant G e ff = n(2 + w)/2. The matter 
density growth factor is primarily sensitive to extensions 
beyond GR in terms of the factor £ = /i(l + zu). These 
parameters still appear to have covariance, however, in 
our initial explorations. Overall, this suggests that ex- 
ploration of gravity through cosmological measurements 
requires a sufficiently flexible theory space and a diverse 
set of observations. 

As seen from Figures [2E1 robust identification of de- 
viations from GR will require measurement over a large 
range of scales. Well below the Hubble scale, the modifi- 
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CJ- (2<z<9) 




FIG. 5: Marginalized probabilities of the post-GR param- 
eters zuoa,b,c defined in high, medium and low redshift bins 
respectively. The parameter fi has been fixed to fx = 1, con- 
sistent with General Relativity. Green (dot-dashed) curves 
are constraints determined from the WMAP 5 year [53] and 
supernova Union2 [6(| data sets only. Red (dashed) curves 
also include the COSMOS weak lensing tomography data [j| . 
Black (solid) curves use measurements of the aperture mass 
taken from the CFHTLS weak lensing survey [4j in addition 
to COSMOS, WMAP5, and Union2. 



cations we have examined become scale independent and 
so can become confused with shifts in the fiducial ampli- 
tude (erg), galaxy bias, or normalization errors from pho- 
tometric redshift estimation of weak lensing source den- 
sities. These will need to be addressed to have confidence 
in claims of any detected deviation, as will allowance for 
expansion histories different from ACDM. 

Finally, future data, including observations sensitive to 
growth and the growth rate, and those sensitive to the 
expansion history, will be essential to providing true tests 




FIG. 6: Marginalized probabilities of the post-GR param- 
eters fioa,b,c defined in high, medium, and low redshift bins 
respectively. The parameter zu has been fixed to zu — 0, con- 
sistent with General Relativity. All curves show constraints 
derived using data from the WMAP 5 year release [5?| . su- 
pernova Union2 set and COSMOS [| plus CFHTLS [| 
weak lensing data. 




FIG. 7: 68% (inner) and 95% (outer) confidence contours 
for modified gravity parameters are plotted in the zuq-^iq 
plane (all other parameters marginalized), where vuo and no 
are defined as in Eq. (|18p . Blue (background) contours use 
the WMAP 5 year, supernova Union2, and COSMOS data 
sets. Red (foreground) contours use WMAP 5 year, super- 
nova Union2, COSMOS and CFHTLS data. The black curve 
plots the degeneracy direction fio = 2/(2 + zuq) — 1 from Fig- 
ure^ The yellow x denotes GR parameter values. 



of the framework of gravity on cosmic scales. 
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Appendix A: Density Perturbation Growth 

We here obtain the analog of the GR second order dif- 
ferential equation for matter density perturbation evolu- 
tion, working in the conformal Newtonian gauge. After 
matter-radiation decoupling, conservation of energy gives 



30 



-ne ri 



(Al) 

(A2) 



assuming 5p = a — 0, i.e. there is no pressure, no pres- 
sure perturbation, and no anisotropic shear. Rearranging 



30 - 5 m 
30 -A r 

30- A r 



d (me. 
rf7 \~~~k 2 
3 r 



k 2 

30 - A m + 3HiJj 
l-3{H~H 2 )/k 2 



H - Hi 



m 2 4> 



(A3) 



We can use Eqs. pip and © to write ip in terms of 
A m , [i, zd, and background quantities; similarly we can 
use the time derivative of Eq. (fTTj) to write 0. This gives 



S m + 



30 



me m me 



me n 

k 2 



m 

k 2 



-ne m + k 2 ^) (A4) 



where the second equality comes from using Eqs. (|A1|) 
and (|X2|) . Substituting Eq. ([M)) into (|A4|) would just 
return the truism A m = A m . However, if we take the 
first conformal time derivative of Eq. (|A4[) before substi- 
tuting, we find a second order differential equation de- 
scribing the evolution of A m for arbitrary w(a, k) and 
/i(a,k). We omit the explicit copious algebra and show 
the result: 



A m I 1 + -2 Tn 



A, 



¥ (2r M + 2Tfi) - ¥ (i + m)v^n + In + 3 fc2 



-3r^ - k 2 



9^ + 3-y X ■ 

k k k 2 -'3[H-H 2 



A m <| - ^ (f/jL + Til + 2f/i) - (1 + m)-£ (-k 2 +6H- 3H 2 ) - Sza^H - ^(1 + zu)U (tfx + T/i) 



< n HH U 3\ -3(r/i + r/i) -3(i + w)wr M ' 

7^ + 3-^-9^^ + 3— I x 



k 2 k 2 



k 2 



k 2 - 3 n - n 
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r 



where T = 4irGa 2 p m . All that we have assumed in this of the background cosmology so that Eqs. (|A1[) and (|A2 
derivation is that matter and A are the only constituents hold. 
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